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Abstract. The properties of polymer liquids on hard and soft substrates are 
^ I investigated by molecular dynamics simulation of a coarse-grained bead-spring model 

Q ' and dynamic single-chain-in-mean-field (SCMF) simulations of a soft, coarse-grained 

O . polymer model. Hard, corrugated substrates are modelled by an FCC Lennard- 

Jones solid while polymer brushes are investigated as a prototypical example of a 
soft, deformable surface. From the molecular simulation we extract the coarse-grained 
parameters that characterise the equilibrium and flow properties of the liquid in contact 
with the substrate: the surface and interface tensions, and the parameters of the 
hydrodynamic boundary condition. The so-determined parameters enter a continuum 
description like the Stokes equation or the lubrication approximation. 
l/^ ' At high temperatures the Navier slip condition provides an appropriate description 

^P , of the flow past hard, corrugated surfaces. The position, xi,, where the hydrodynamic 

boundary condition is to be enforced, agrees with the location of the liquid-solid 

interface and the slip length can be consistently identified by comparing planar shear 

flow and parabolic, pressure-driven flow. If the surface become strongly attractive or 

^^ ' the surface is coated with a brush, the Navier slip condition will fail to consistently 

describe the flow at the boundary. This failure can be traced back to a boundary layer 
with an effective, higher viscosity. 

The solvent flow past a polymer brush induces a cyclic, tumbling motion of the 
tethered chain molecules. The collective motion gives rise to an inversion of the flow in 
the vicinity of the grafting surfaces and leads to strong, non-Gaussian fluctuations of 
the molecular orientations in the flow. Both, molecular dynamics as well as dynamic 
SCMF simulations, provide evidence that the flow past a polymer brush cannot be 
described by Brinkmann's equation. 

The hydrodynamic boundary condition is an important parameter for predicting 
the motion of polymer droplets on a surface under the influence of an external force. 
The steady state velocity is dictated by a balance between the power that is provided by 
the external force and the dissipation. If there is slippage at the liquid-solid interface, 
the friction at the solid-liquid interface and the viscous dissipation of the flow inside 
the drop will be the dominant dissipation mechanisms; dissipation at the three-phase 
contact line appears to be less important on a hard surface. 

On a soft, deformable substrate like a polymer brush, we observe a lifting up of 
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the three-phase contact Hnc. ControlHng the grafting density and the incompatibihty 
between the brush and the polymer hquid we can independently tune the softness of 
the surface and the contact angle and thereby identify the parameters to maximise the 
deformation at the three-phase contact. 
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Figure 1. W||(a;), tangential to the solid substrate. Note that the velocity of the 
liquid never reaches zero, i.e., there is an apparent, microscopic slip. The red dashed 
lines indicate the linear and parabolic hydrodynamic velocity profiles, i;,, ' -^ '^° and 

1^11 ' ■*' ^° , according to Eq. (^5)) and Eq. (PH) for Couette and Poiseuille flow, respectively. 
These profiles are obtained by fitting the measured velocity profiles, vAx)^ far away 
from the solid substrate, to the asymptotic hydrodynamic continuum prediction. The 
extrapolation of these hydrodynamic velocity profiles to zero mark the positions, xq, 
and xp. The green triangle indicates the Navier-slip condition, Eq. ([2]), and the 
hydrodynamic position, Xb, and slip length, (5, are indicated. 



1. Introduction 

The equilibrium properties of a liquid on a solid substrate are chiefly determined by 
the balance of surface and interface tensions [H El El HI |5l E]. Balancing the tensions 
parallel to the substrate at the three-phase contact between the solid, the liquid and its 
coexisting vapour, Young derived in 1805 the relation [1] 

7lv cos 6 + 7ls = 7vs (1) 

between the macroscopic contact angle, 6, of a drop, the liquid-vapour interface tension, 
7lv, and the surface tensions, 7ls and 7vs5 of the liquid and its vapour in contact with 
the substrate. These tensions depend on the microscopic details of the liquid and the 
substrate, and their prediction requires molecular simulation of a microscopic model 
of the materials. Using the so-determined tensions, one can use Young's equation to 
predict the macroscopic equilibrium properties of droplets on surfaces. 

The flow of liquids past surfaces is of paramount importance in engineering 
applications [3 El 13 |TOl |TT]. On a macroscopic scale it is often assumed that the 
liquid velocity at a boundary equals the velocity of the surface. On a microscopic scale, 
however, this no-slip condition can be violated and the molecules of the liquid may 
slip past the surface. This slippage at the solid-liquid interface and the concomitant 
details of the velocity profile at the solid-liquid boundary have attracted abiding interest. 
Balancing the viscous stress due to the velocity gradient inside the fluid with shear 
viscosity, r^, against the friction stress of the liquid slipping past the solid substrate. 
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Navier formulated in 1823 a hydrodynamic boundary condition [12 

dv 



Xvn 



Xb dx 



(2) 



friction stress at solid 



Xb 

viscous stress in fluid 



where A is a friction coefficient that quantifies the hnear relation between the tangential 
velocity at the surface v\\\x^^ and the friction stress. Xb denotes the position, at which 
this Navier-slip boundary condition is to be applied. The ratio 6 = rj/X is a length 
and is denoted as slip length as illustrated in figure [U This condition serves as a 
boundary condition to a continuum description of the fiuid fiow by a Navier-Stokes 
equation. Molecular simulations can measure the viscosity, rj, of the bulk liquid and 
the two parameters - the shp length, 6, and the position, Xb, of the hydrodynamic 
boundary condition - of the Navier-slip condition and thereby transfer information from 
a microscopic, particle-based description to a continuum one. 

In the following we will utilise molecular simulations of coarse-grained polymer 
models to extract the parameters that relate the fiow properties of a particle-based 
model to a macroscopic description and test the validity of the assumptions of the 
macroscopic description on small length scales. In the next section we will introduce 
the particle-based models and simulation techniques. Then, we investigate the fiow 
past hard, corrugated substrates and soft, penetrable, brush-coated substrates. We find 
that the routinely used Navier-slip boundary condition may fail to describe the near- 
surface velocity profiles independent from the type of fiow [131 [IH US] cind rationalise 
this observation by a simple, phenomenological model [H]. Moreover, we observe the 
reversal of the fiow direction in the vicinity of the grafting surface of a polymer brush 
in fiow [13]. The fourth section studies the behaviour of polymer drops on hard and 
soft surfaces and illustrates that the slip at the solid-liquid boundary is an important 
dissipation mechanism for the motion of small drops [16]. The manuscript concludes 
with a brief outlook on dewetting of thin liquid films and perspectives to future molecular 
simulations. 

2. Models and Methods 

The surface and interface tensions as well as the viscosity and friction of the motion of 
simple liquids at surfaces ultimately depend on the molecular structure on the atomistic 
scale. In polymer systems, however, the extended structure of the molecules imparts 
some universal aspects onto the behaviour of these macromolecules at surfaces and 
interfaces. Varying the molecular weight and, in case of brush-coated surfaces, the 
grafting density one can experimentally control the wettability, adhesion, and friction 
of a surface and the shear viscosity without altering the interactions on the atomistic 
scale. Therefore, polymeric systems are often described by coarse-grained models. In 
the following, we will use two distinct computational approaches - molecular dynamics 
simulation of a Lennard- Jones, bead-spring model [TTJ [181 [IS] ci-nd dynamic single-chain- 
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in-mean-field (SCMF) simulations of a soft, coarse-grained model of a dense polymer 
melt [201 1211122]. 

2.1. Molecular dynamics of a Lennard- Jones bead-spring model 

Molecular dynamics simulations are performed using a widely used coarse-grained 
polymer model [U]. A macromolecule is treated as a string of A^ beads, connected by 
nonlinear springs to form a linear chain molecule. The beads interact with a Lennard- 
Jones potential 

f/„tW=4.„4(^)"-(^)'] (3) 

For particle distances, r, larger than a cut-off distance, re, the potential is cut-off and 
shifted such that it is continuous at r^. r^ = 2^/®o"q,^ gives rise to a purely repulsive 
interaction [17J, which is appropriate to simulate dense polymer liquids. In order to 
study polymer droplets with a free surface, the liquid must coexist with its vapour 
and, to this end, the non-bonded interactions must include a longer-ranged, attractive 
contribution. This attraction can be included by choosing re = 2 ■ 2^/^(Tq/3 [iHl [El [23] or 
Tc = 2.5aai3 [21]. Indices, a,/3, are related to different bead species, i.e., segments of the 
polymer liquid, segments of the polymer brush, or particles of the substrate, e and a 
set the energy and length scales of the Lennard- Jones model. Adjacent segments along 
the chain molecule are connected by an anharmonic, finite extensible, nonlinear elastic 
(FENE) potential [IT], 

2" 



U^ir) = -hRl\n 



. r 
Ro 



(4) 



with k = 30e/(7^ and Ro = l.5a. These parameters are chosen so that unphysical bond 
crossings and chain breaking are eliminated. In the following, all quantities will be 
expressed in terms of molecular diameter a = 1, binding energy e = 1 and characteristic 
time r = y^ma'^/e. 

The equations of motion are solved via the velocity- Verlet algorithm. The 
temperature is kept constant using a dissipative particle dynamics (DPD) thermostat 
[251 [26l [271 [28l [29] . This thermostat conserves the local momentum and thus results in 
hydrodynamic flow behaviour on large length and long time scales. It adds to the total 
conservative force, Eq. ([3]) and Eq. (j4]), a dissipative force, F^ , and a random force, 
F^ . Both forces are applied in a pairwise manner, such that the sum of thermostatting 
forces acting on a particle pair vanishes. Let F be the friction constant, the dissipative 
force takes the form 

Ff = -r5^a;^(r,,)(r,,.v,,)r,, (5) 

where Vij = (rj — Yj)/rij and Vjj = Vj — Vj. We choose the weight functions: 

a;^(r,,) = {^^~ ''"'^'''^^ '''' ^ "'" (6) 

^ " 1 r^,- > r, ^ ' 
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N 


10 


100 


200 


-Reo [O-] 


3.57 


12.28 


15.2 


T] [y/me/cr^\ 


6.4 ±0.1 


73 ±4 


120 ±6 


6D [aVr]xl03 


78 


3.1 


1.6 



Table 1. 

PLCT^ « 0.8 



Bulk properties of the Lennard- Jones polymer liquid for fc^T = 1.2e, 
and rr — 2.5a. 



with Tc identical to the one used in Eq. (|3]). The random force is given by: 



riR 






u 



R, 



rij)0, 



ijT^ij 



(7) 



where 6ij is a random variable with zero mean, unit variance, second moment 



{9ij{t)dki{t')) = {6,,6,i + 6ii6jk)6{t-t'), and % = 
satisfy the fluctuation-dissipation theorem, [w^] 



6.ji. The weight functions, u 



Ri 



^jji 



u 



D 



Friction, F, and noise strength 



^ define the temperature via ^^ = 2kBTT. We choose F = 0.5/r in all our simulations. 
The equation of motions are integrated with a time step of At/r = 0.002 or 0.005. The 
simulations are performed using a parallel molecular dynamics program based on force- 
decomposition [30], which is particularly suitable for spatially inhomogeneous systems 
like droplets, or LAMMPS [31], which employs a geometric parallelisation strategy. 

Some bulk properties for cut-off r^. = 2.5a are compiled in tabled] The data for the 
model with re = 2-2^/^cr are similar. Two temperatures, ksT/e = 1.68 [32ll33l[34l[35l[T5] 
and 1.2 [161 [E]) have been studied for r^ = 2 • 2^/^cr and the density of the liquid that 
coexists with its vapour is plO"'^ = 0.61 and 0.79, respectively. The density of the vapour 
is neghgible at these temperatures far below the Theta-temperature, fcsG/e ~ 3.3, of our 
model [19]. The spatial extension of a chain molecule is characterised by an average, 
mean-squared, end-to-end polymer distance R^^ = (i?g)^/^. For A^ = 10 we observe 
-Reo/c = 3.7(1) and 3.6(1) at the higher and lower temperatures, respectively. Thus, this 
chain length corresponds to very small values of the invariant degree of polymerisation, 
A/" = {pi^RlJNf = 9 and 14. 

The surface tension, 7, of the liquid-vapour interface can be measured using a slab 
geometry [37] [191 EH] where the interfaces are parallel to the xy-plane. 7 is determined 
by the anisotropy of the pressure 

(P..) + (Pyy) 



2^A = V[ (P,,) - 



(8) 



where the interfacial area is A = L^ x Ly and the factor 2 arises because the periodic 
boundary conditions stabilise two liquid-vapour interfaces in the simulation box. The 
components of the pressure tensor are given by virial formula 



a/3 



\ i i<j / 



(9) 
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Figure 2. ksT/e = 1.2 maintaining a steady-state momentum flux by exchanging of 
particle momenta in the shaded region according to Ref. 36i . Ns denotes the number of 
molecular dynamics steps between the exchanges of momenta. The inset extrapolates 
the results to vanishing momentum flux. The cubic simulation cell has a spatial extent 
of L = 29.3(7 and contained 19 160 Lennard- Jones segments. 



where the sum over i and i < j runs over all particles or pairs of particles, respectively. 
m is the mass of a segment, v its velocity and Fjjq denotes the component, a, of the 
force between particles i and j that are a distance rij apart. For ksT /e = 1.2 this 
procedure yields 7cr^/e = 0.515(3) [16j, while grand-canonical Monte-Carlo simulations 
yield the value ja'^/e = 0.16(1) for fc^T/e = 1.68 and r^ = 2 ■ 2^/^a [19]. 

The Rouse time tr is related to the self-diffusion coefficient via D = R^^/ (Svr^rp). 
For the high temperature one finds D = 0.05o"^/r (for Tc = 2 ■ 2^/®cr), and at the lower 
temperature one obtains D = 0.0157(30)o-Vr or 0.013(3)0-7^ for r^ = 2 ■ 2^/^a and 
Tc = 2.5a, respectively. The Rouse time characterises the longest relaxation time of 
the molecular conformations in equilibrium. The product of tr and the shear rate, 7, 
defines the Weissenberg number. 

The shear viscosity, rj, has been computed using the reverse, non-equilibrium 
molecular dynamics method |36t |39[ HO]. To this end, one divides the system into 
slabs and exchanges the momenta between the two slabs that are separated by half the 
system size as indicated by the shaded regions in figure [2l At each momentum swap, 
one exchanges the tangential momenta, Py = mv\\ , of the particles with the smallest v\\ 
in the centre, shaded slab with the momenta of the particle with the largest v\\ in the 
left, shaded slab. This procedure imposes a momentum transfer between the slabs and 
results in a linear velocity profile, v\\{x), between the shaded regions. From the gradient 
of the velocity and the momentum fiux, one obtains the viscosity according to [36] 



V 



^xy 



dx 



with S 



xy 



Ap/At 
2A 



(10) 



where Ap denotes the momenta exchange during the time interval. At, and J^^y the 
shear stress. The number of exchange momentum swaps Ng = 3 and Ng = 120 were 
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chosen in order to verify a linear velocity profile inside the fluid. Depending on A^^, 
we performed simulation runs of length between 5 OOOr and 24 OOOr, while data were 
collected after half the simulation run, when a steady state is reached. The simulation 
cell was comprised of 19 160 particles. These results for the viscosity of the model with 
cut-off Tc = 2.5(7 are very similar to the results rj = 5.3(2) ^Jrakja^ p^ obtained from a 
Green-Kubo relation 



^GK(t) = -^jdt' (P.,(t')P.,(0)) (11) 

for the model with Vc = 2-2^/^0" at ksT /t = 1.2 using a molecular dynamics run of length 
200 OOOr for 5440 particles [16j. At the higher temperature, ksT/e = 1.68, measuring 
the stress in a non-equilibrium molecular dynamics simulation yields t] = 1.9(3) y/me/a^ 
[Mf [T5] . In the Lennard- Jones bead-spring model the friction coefficient of the DPD 
thermostat is chosen small enough that its value does not have a pronounced influence on 
the dynamics, i.e., the segmental friction chiefly stems from the bonded and non-bonded 
interactions. 

2.2. Dynamic SCMF simulations of a soft, coarse-grained polymer model 

Additionally to the standard, bead-spring model we also employed a soft, coarse-grained 
model of a dense polymer liquid [20l [2T1, [22] . In this model, neighbouring segments along 
the backbone of a chain molecules are bonded together via simple, harmonic springs of 
the form 

where i?eo denotes the average end-to-end distance. 

Non-bonded interactions give rise to a small compressibility of the dense polymer 
liquid and are described by [42] 

^nb = ^PL/d^r (0^(r)-l)^ (13) 

where 0(r) denotes the normalised density of the liquid which is calculated from the 
explicit particle coordinates. We choose the value KqN = 50, which is sufficiently large 
to suppress density fluctuations on the length scale of polymer coils. These density are 
efficiently calculated via a collocation grid. The quadratic form of Tinh can be rewritten 
as a sum over pairwise interactions, f/nb, between segments [21]. In contrast to the 
Lennard- Jones potential, however, these pairwise interactions do not exhibit a harsh 
repulsion at short distances but rather correspond to a soft, repulsive interactions that 
allows coarse-grained segments to overlap. Qualitatively, this soft, repulsive interaction 
between effective segments, which each represent a collection of atoms, would be the 
expected outcome if the interactions were derived by explicitly integrating out the atomic 
degrees of freedom (systematic coarse-graining) [lU 
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In a Lennard- Jones, bead-spring model the monomer density cannot significantly 
increase above pi^a^ ~ pi^R^^/N^/"^ ~ 0{1) because the fluid of segments either 
crystallises or becomes a glass at higher densities. In a soft, coarse-grained model, 
the segment density can be increased much further, pi^R^^/N^/"^ ~ 10 — 100, and 
thus the model allows us to assess chemically realistic values of the invariant degree 
of polymerisation, M = [pi^Rl^/N]^ = [p-^Iil^/N'^/'^YN can be studied by molecular 
simulation. In section. 13.31 we will study molecules with AT = 16 384, which is not 
feasible with a Lennard- Jones bead-spring model [45] . 

The properties of the soft, coarse-grained model are studied by dynamic SCMF 
simulations. In SCMF simulations [211 HHl E] one considers a large ensemble of explicit 
chain configurations that independently evolve in an external field. This field 

W{r) = -^ (14) 

PL d(j){r) 

mimics the effect of the non-bonded interactions. Its gradient gives rise to a force 

In dynamic SCMF simulations, chain segments evolve via self-consistent Brownian 
dynamics [Ml |49l [50] or Smart-Monte-Carlo (SMC) moves [SH |52] with trial 
displacements 



Artnai = (V) At + ^At + eW ^^^ (15) 

Co y s>o 

(v(r)) denotes the hydrodynamic velocity field. F = Fb + Fnb + Fgx is the force acting on 
a segment arising from bonded and non-bonded interactions as well as external potentials 
due to confining surfaces, (o characterises the segmental friction, and ^ is a Gaussian 
random number with zero mean and unit variance. Such a trial displacement, r — ;■ r, is 
accepted with probability 



A{r —7- r') = min I 1, exp 



1 



ksT 



E{r') - E{r) 

(16) 



F(r') + F(r),, , AA , , „, ^, ,o. 
- \ W -^} + -r {F(rT - F(r)2} 



where E = Xlbonds Uh + 'Hnh denotes the total energy of the system. The SMC-algorithm 
samples the equilibrium properties independent from the choice of the forces, F, used 
to bias the proposed trial moves or the magnitude of At. In the limit At — )■ 0, 
the acceptance rate of the SMC moves tends to unity and the algorithm reproduces 
Brownian dynamics. One can use however much larger time steps. At = 0.08 ^^^ °^ 
than in Brownian dynamics and still reproduce the over-damped, diffusive motion 
of the segments in a dense fluid ^3]. In this case, the self-diffusion coefficient and 
the shear viscosity in the bulk are close to the Rouse limit, Z)r = k bT /{CoN) and 
^R-Reo/(Co-^v-^) = 1/36, respectively. The acceptance rate of SMC moves is about 
83%. Chains exhibit Rouse-like dynamics, which is characteristic for unentangled melts 
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Figure 3. L^ — lOReo and the lateral dimensions are Ly = L^ = 4_Ro 
normalised density, </>, is comprised of the density, <pai, of the melt and the density, 0b, 
of the grafted chain molecules. Data for Poiseuille flow (/|| A^ = O.OSfcsT/i^eo, symbols) 
and Couette flow {jo^oNR'^^/kBT — 5, Wio — 0.247, lines) are shown. The simulation 
cell contains 655 360 segments, b) Kramer's stress for Poiseuille and Couette flow. 
The horizontal dashed line marks the stress obtained from the forces on the grafted 
ends for Couette flow. Adapted from Ref. [I3] 



and the (bare) Weissenberg number is defined as Wio = %CoN R'^^/^Sn'^kBT) where 
7o denotes the shear rate imposed by moving the confining surfaces with the grafted 
polymer chains. 

In the Lennard- Jones, bead-spring model, the non-bonded interactions provide an 
important contribution to the friction and viscosity of the liquid of monomers gives rise 
to corrections to the chain-length dependence of the viscosity predicted by the Rouse 
or reptation model. In the soft, coarse-grained model, however, the the non-bonded 
interactions hardly give rise to friction. Instead, the friction is chiefiy generated by 
the random forces in case of Brownian dynamics, i.e., the parameter, (o, of the SMC 
algorithm dictates the friction. Moreover, the soft non-bonded potentials in conjunction 
with the harmonic bonded ones allow the chains to cross through each other in the 
course of their motion in SCMF simulations [53]. 

The fluctuating, external field, W, is calculated after every time step. It is important 
that they are frequently updated such that they mimic the instantaneous, non-bonded 
interactions. If this condition is fulfilled [21] SCMF simulations will capture even subtle 
deviations from the Gaussian statistics in a dense polymer melt [55] • The hydrodynamic 
velocity, (v), however, represents the average flow field and must not fiuctuate. It is 
calculated self-consistently from the particle displacements. First, we directly calculate 
an instantaneous velocity, v, = Arj/ At, of segment i from its explicit displacement [56] 
during a SMC step (in order to retain spatial resolution) and assign it to the grid. Then, 
we add the instantaneous velocities of all segments, average over a time period, T, and 
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normalise by the local density to obtain the average local velocity (v(r)). This procedure 
ensures that the average force, (F), vanishes. This time averaging procedure limits the 
simulation technique to stationary or slowly varying flows. For the density utilised in our 
simulations a time interval T > 800(oRlo/i-^^BT) is sufficient to eliminate fluctuations 
of the velocity field and to yield accurate predictions in equihbrium. 

In figure [3^ we present the density profile across the slit pore with width, L^ = 10i?e- 
The solid surfaces are coated with a polymer brush and the tethered chains have the 
same physical properties as the chains in the melt. The grafting density is (JgRl^ = vAf. 
Results of dynamic SCMF simulations are presented for Couette and Poiseuille flow. At 
this moderate grafting density, there is a broad interface between brush and melt ("wet 
brush"). In agreement with previous self-consistent Brownian dynamics studies j49] 
and molecular dynamics simulation of bead-spring models [241 El] , the profile normal 
to the surface does not depend on flow for small shear rates. The depletion at the 
wall stems from a repulsive segment-wall potential U^a.n{x) = A^ exp{—x'^ /2e^^i{) with 
A„A^ = 66.6^^7" and ew = 0.15i?e- This repulsive potential impedes the Smart-Monte- 
Carlo algorithm to propose moves that would penetrate the hard walls and be rejected. 

Measuring the forces, F.f'^^ , that act on the grafted segments in Couette flow, we 
obtain the shear stress 

^igraft p3 

0.229VAr (17) 



ksT LyL^ksT ksT 



for Couette flow. This value is indicated by the horizontal dashed line in figure lib) and 
it agrees well with the prediction of the Rouse model 

yRouse p3 2 

"^ "^ - Wio ^ 0.203 (18) 



fcsTvW 12 

Alternatively, we can estimate the intramolecular stress of the melt at the centre 
of the film via the mean- field approximation (Kramer's formula [54]z.e., virial of the 
bonded interactions, Eq. f fT2|) ). which only accounts for the bonded contribution to the 
stress and disregards the contribution of all non-bonded forces. 

yKrainer p3 Ih h \ 

-^ -^ = [N -\U^--ppil (19) 

where hr^ and hy denote the distance between bonded segments perpendicular to the 
surfaces and along the shear direction, respectively. This estimate is shown in panel (b) 
of figure [3] for Poiseuille and Couette flow. We have utilised the same spatial assignment 
for the stress as for the density. Alternative assignment schemes for the stress that are 
locally more accurate can be envisioned [57] . Thus the small variation of the stress in 
Couette flow may be either due to inaccuracies of the local assignment or contributions 
that stem from the non-bonded forces suppression density fluctuations. 

The stress in Couette geometry agrees well with the result obtained from the 
force on the grafted segments. From the stress at the centre of the film and the 
velocity gradient (see figure EH) we can estimate the dimensionless shear viscosity. 
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Figure 4. es = 0.6e = 0.73eg°* and temperature, ksT/e — 1.2. The measured velocity 
profiles are marked by symbols, the hydrodynamic velocity profiles obtained by fitting 
the simulation data at the centre of the film at indicated by red, dashed lines and 
the construction of the slip length, S, and the hydrodynamic position is depicted by 
the dashed-dotted, green lines. Inset: Slip length, 5, as a function of the attraction, 
0.2 < es/e < 1, between solid and liquid at temperature, fc^T/e = 1.2. The solid line 
with circles is obtained from the Couette and Poiseuille profiles using non-equilibrium 
molecular dynamics (NEMD) simulation, while the dashed line with squares from the 
Green-Kubo (GK) relation Eq. ((20|) . The dashed line marks the behaviour S ^ l/cg 
according to Eq. (|2ip . Adapted from Ref . [16] . 



SGrjUeo/iCoNVAf) = 1.12 for Couette flow. Tlie linear dependence of the stress in 
Poiseuille flow shown in figure [3]d is well described by —^^ = O.OAlAkBTvW / R^^ and 

the curvature of the parabolic velocity profile (cf. figure[7K) is -^^ = 1.38kBT/{(oNRl 
Both results yield 36 r]R^o/{CoNV^) = 1.08. 



^e^o) 



3. Flow past hard and soft substrates 

3.1. Corrugated, hard substrates 

First, we consider the flow of a Lennard- Jones, bead-spring liquid past a hard, corrugated 
substrate [I6], which is comprised of two layers of an FCC lattice with density psC"^ = 2. 
The particles of the solid substrate interact with the segments of the liquid via a Lennard- 
Jones potential with as = 0.75 a and es = 0.2 — 1.0 e and Tc = 2 ■ 2^/^o". Increasing es we 
simultaneously alter the wettability of the surface in equilibrium as well as the friction. 
The friction coefficient, A in Eq. ([2]) can be evaluated via a Green-Kubo relation 

mmM 



A 



1 



dt (FLs|lWi^LS||(0)> 



(20) 



ksTA Jo 
for the tangential forces, -^LSlh that the substrate exerts onto the fluid. A denotes 
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the surface area. Using the shear viscosity, rj, we plot the shp length S = rj/X in the 
inset of figure H] for ksT /e = 1.2 and various values of the attraction, es, between 
substrate and liquid. The slip length decreases as we increase the attraction between 
liquid and substrate, but 6 remains finite at the wetting transition, which is located 
around €3^* ~ 0.82e. Qualitatively, this finding agrees with other simulation studies 

[601 iH [621 ESI El- 

Barrat and Bocquet devised an approximation [651 [H] 



which captures the main ingredients that qualitatively dictate the friction at the surface. 
S{q\\) is the structure factor of the liquid at the surface and q\\ is the wave vector of the 
surface corrugation. -D(g|j) denotes the collective diffusion coefficient at wave vector, gy. 
Assuming that the normalised density profile (j){z) is not altered by the strength es of 
the interaction, [/ls, between liquid and substrate, one obtains the prediction, 6 ~ l/cg, 
which is nicely confirmed by the simulation data in figure HI Moreover, the adhesion free 
energy W = 7lv(1 + cos 6) roughly depends linearly on es for es > 0.4e [16], because 

des A es 

with Tis being the energy of the interaction between the sohd substrate and the hquid, 
is chiefiy dictated by the normalised density profile, 0(x). The latter, however, only 
weakly depends on es provided that (i) the temperature is sufficiently low and the liquid 
compressibility is concomitantly low and (ii) es is sufficiently large in order to suppress 
the formation of a drying layer. In this regime, the data are also compatible with a 
behaviour of the form, 5 ~ (1 + cos0)~^ , which has recently been observed for water 
at hydrophobic surfaces [66l |6^ . 

Alternatively, one can determine the two parameters - 6 and Xb - of the Navier 
slip condition by non-equilibrium molecular dynamics. Using Couette fiow, one can 
only extract the combination, Xb + S. Note that for chemically complex surfaces or 
soft, deformable substrates {e.g., brush-coated surfaces or networks), the hydrodynamic 
position, Xb, of the substrate in the continuum description may not be obvious, and it 
may depend on external control parameters (like temperature, pressure and grafting 
density) [341 [15]. In order to determine both parameters independently, we study 
Couette and Poiseuille fiow. These two fiow profiles are particularly simple in the sense 
that they are the only ones, for which the stress depends only linearly on the distance 
from the surfaces [68]. Therefore, a non-locality of the viscosity, which is often invoked 
by generalised hydrodynamic descriptions [691 HO] in order to extend the continuum 
description to smaller length scales, has no pronounced consequences. 

Couette fiow in a thin film between two parallel, hard, corrugated walls is produced 
in non-equilibrium molecular dynamics simulations by moving the confining surfaces in 
opposite directions with velocity v\\ = ±70/ (2L^), where L^ denotes the film thickness. 
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This results in a linear velocity profile v., ' ^ ^°{x) at the centre of the film, where a 
hydrodynamic continuum description is valid. 

^C,hydro^^^^^^^_^^^ (23) 

At the surface, however, there may be deviations from the linear velocity profile as 
illustrated in figure [1] From the velocity profile at the centre one extracts the shear 
rate, 7, which may deviate from the value % imposed by the moving surfaces and the 
position, xq, where the hydrodynamic velocity profile extrapolates to zero. 

Poiseuille fiow in a thin film is generated by applying a force, f\\, onto all segments. 
This results in a parabolic velocity profile, v,, ' ^ ^°(x), at the centre of the film, which 
takes the form 

P,hydro^^^ = M (a; - Xp) (L. - Xp - x) (24) 

Zt] 

Provided that the film thickness, L^, is large enough to observe bulk-like behaviour at 
the centre, the non-equilibrium simulation of Poiseuille fiow yields the shear viscosity, rj, 
in the bulk and the position, xp, where the hydrodynamic velocity profile extrapolates 
to zero. These extrapolations are presented in figure H] for es = 0.6e. 

Since the Navier-slip condition provides a boundary condition to the hydrodynamic 
continuum description one should use the extrapolations of the velocity profiles at the 
centre to extract the slip length, 6, and hydrodynamic position, x^. These parameters 
are simply related to xc and xp via [HI [M] 



6 = a/(xc - xp){L^ -Xp - Xc) (25) 

2:b = a^c + (^ (26) 

The results of these non-equilibrium molecular dynamics simulations are also 
plotted in the inset of figure HI Gratifyingly they agree with the data extracted from the 
Green-Kubo relation for small and moderate attraction between polymer and substrate. 
For strong attraction, e^, the Green-Kubo method however yields slightly larger values 
of the slip length than the non-equilibrium method. This deviation is partially due to 
the fact that we employ the bulk value of the shear viscosity while the structure and 
dynamics of the liquid is altered by the strong interactions with the substrate [71] . 

The temperature dependence of the slip length, 6, and the hydrodynamic position, 
Xb, extracted from non-equilibrium molecular dynamics simulations are presented in 
figureO Data sets correspond to various strengths, es, of the attraction between the solid 
and liquid as indicated in the key, and the liquid density corresponds to the temperature- 
dependent coexistence value. If the surface is only weakly attractive, es = 0.4e, the slip 
length, 6, and the hydrodynamic position, Xb, exhibit only a very weak temperature 
dependence. At low temperatures, however, when we approach the glass transition of 
the model [72l[73] at around fc^Tgiass/e ~ 0.41, the slip length, 6, increases. The growth 
oi 6 = rj/X for T — )■ Tgiass indicates that the friction inside the liquid, rj, diverges more 
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Figure 5. S, versus temperature for different strengths, es, of solid-fluid interaction. 
The inset represents the position of the hydrodynamic boundary, Xb, versus 
temperature. All lengths, S and Xh, are measured in units of the Lennard- Jones 
parameter, ctlj, and temperature, T, in units of the glass transition temperature, 
fcsTgiass — 0.41e. The thickness of the film is varied between 8 < D/Rco < 17, the 
typical velocity of the wall in Couette flow is Vg = 0.2a /t and the typical volume force 
applied to generate Poiseuille flow is 0.001 < /||c/e < 0.005. Adapted from Ref. [16] 



rapidly than the friction, rj, of the hquid with the corrugated substrate, whose structure 
is independent of T. For T < Tgiass, the polymer systems will become a glass sohd and 
will move like a solid body (plug flow), which corresponds to an infinite value of the slip 
length. Also the hydrodynamic position, x^, increases slightly for T — )■ Tgiass due to the 
increase of liquid structuring (packing) at the substrate upon cooling the system. 

Already at ksT/e = 0.5, i.e., about 20% above the glass transition temperature of 
our model, 6 has increased by an order of magnitude compared to the approximately 
constant value at high temperature. This observation offers an explanation for the 
surprisingly large slip length observed in the dewetting experiments of Fetzer et al. 
[74l [75l [76l [77] , which were performed in the vicinity of the glass transition temperature. 

If the attraction between solid and liquid becomes larger, we observe that the slip 
length, 5, will first decrease upon cooling, pass through a minimum, and eventually it 
will diverge upon approaching Tgiass from above for the reasons explained above. The 
decrease of 5 upon cooling from large T can be rationalised by the changes of the 
liquid structure in the vicinity of the solid substrate. Upon cooling, the liquid becomes 
more structured and this gives rise to an increased, effective viscosity in the vicinity 
of the surface [H]. Quahtatively, decreasing the temperature at constant es is similar 
to increasing the attraction, es, at constant temperature; both effects lead to a more 
pronounced structuring of the liquid at the surface and a smaller slip length. 

At large attraction between solid and liquid, we find that x-p > xq and Eq. ( l25ll has 
no solution. The individual velocity profiles observed for Couette and Poiseuille flow in 
the non-equilibrium simulations do not exhibit any unusual features compared to the 
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Figure 6. Ss — and the formation of a lubrication layer with 77s < 77 at the 
solid substrate. The velocity profiles for Poiseuille (circles) and Couette flow (black, 
solid line), the hydrodynamic velocity profiles at the centre and its extrapolation (red, 
dashed lines), and the construction of the slip length, d, and hydrodynamic position, 
Xb, (green, dashed-dotted line) are sketched. 



situation, where a slip length can be extracted. This observation marks the failure of 
the Navier-slip condition, Eq. (J2]), to consistently describe both types of flows with the 
same set of parameters, 6 and x^. In this situation, Eq. (|2]) cannot be utilised as a 
boundary condition for a continuum description because a boundary condition should 
parameterise the microscopic phenomena at the boundary independently from the type 
and strength of the hydrodynamic flow in the bulk. 



3.2. Failure of the Navier-slip condition: a schematic, two-layer model 

Couette and Poiseuille flows are characterised by constant or linear stress profiles, 
respectively. Thus, we can rule out that the failure of the Navier-slip condition stems 
from a non-local relationship between velocity gradient and stress that might arise from 
the extended molecular conformations of the chain molecules. Nevertheless, the non- 
Newtonian nature of the polymer liquid may contribute to the failure because a finite 
driving must be applied in the non-equilibrium simulations. Pronounced effects of shear- 
thinning on the slip length in polymer fluids have been observed [63] IM]. Moreover, due 
to the spatial extension of the chain molecules, limited range of film thicknesses, L^, 
accessible in the simulation may not be sufficient to observe bulk-like behaviour at the 
centre of the film. 

Therefore, it is beneficial to rationalise the failure of the Navier-slip boundary 
condition in a schematic model, which qualitatively reproduces the failure observed 
in the non-equilibrium simulations but does not invoke the additional molecular length 
scales, a and Reo, strictly deals with a Newtonian fiuid, and can be analytically studied 
for all film thicknesses, L^. In order to capture the changes of the liquid in the vicinity of 
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the solid substrate, we consider a two-layer model, which is illustrated in figure [6l where 
we crudely approximate the gradual variation of the fiuid properties as a function of 
the distance, x, from the solid surface by a boundary region of width, A. For distances 
larger than A the fiuid posses the bulk properties and it is characterised by a viscosity, 
rj. Closer to the solid substrate, a; < A, the hquid properties deviate from the bulk 
and, for simplicity, we parameterise these changes by a constant effective viscosity, ?7s. 
The width is the boundary layer. A, is the only length scale in the model. Within the 
boundary layer, < x < A, and in the bulk, x > A, the Newtonian liquid is described 
by the Navier-Stokes equation. The tangential velocities are denoted by vs{x) and v{x) 
in the boundary layer and the bulk, respectively. At the interface between the solid 
surface and the boundary layer, x = 0, we impose a Navier-slip condition, Eq. (|2]), with 
a microscopic slip length, ^s- At the interface between the boundary layer and bulk, we 
require the continuity of shear stress and velocity. 



dvs 
dx 



dv 
dx 



and ^^sL=A- = ^L=A+ (27) 

=A+ 



These are the natural boundary conditions at a fluid-fluid interface. 

One can straightforwardly calculate the flow proflles in the boundary region and in 
the bulk. For Couette flow, one obtains: 



Tj 

fs(x) = 7 — ix — (5s) for x < A (2^ 

^s 



f (x) = 7 (x — xc) for X > A (29) 



with Xc = ( 1 - — I A - —5s 



A similar calculation yields 



vs{x) = - ^ (x^ - L^x - D5s) for x < A (30) 



2r/s 

h 

2ri 



v{x) = (x — xp) (L^. — xp — x) for x > A (31) 



with Xp (Lx — Xp) = — I 1 1 A^ + LxXc 

\ VsJ 

for Poiseuille flow. L^. denotes the fllm thickness, which is deflned by the distance 
between the positions where the microscopic Navier slip condition with slip length, 6s, 
is applied. 

Finally, Eq. ( l25l) yields the slip length 



A^(^-l]{A + 26s)+(^6sy (32) 

which is independent from the fllm thickness, L^.. In this schematic model, the order of 
magnitude of the slip length is set by the spatial extent. A, of the boundary layer. 
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The first term describes the effect of the surface layer, the second term arises from 
the microscopic slip at the solid surface. The condition, rj/rjs > 1, corresponds to the 
formation of a lubrication layer at the surface and results in an enhanced slip length, 
S > Ss, compared to the microscopic slip at the solid- fluid interface. If, however, the 
solid-liquid interactions give rise to a boundary layer with a large effective viscosity, 
rj/rjs < 1) the presence of this sticky boundary layer reduces the slip length, 6 < 6s- 
Moreover, if 

, < 1±2M^ (33) 

the Navier-slip condition, Eq. ([2]), fails and xp > xq as observed in the simulations of a 
polymer liquid in contact with a strongly attractive substrate. 

Qualitatively, our schematic model can rationalise the observations in our molecular 
simulation: (i) At high temperature, kinetic effects will dominate the behaviour, thus 
rjs ~ rj. In this case, 6 is equal to the microscopic slip length 6 ^ 6s- (ii) Upon cooling 
the fluid, the shear viscosity increases. If the solid-fluid interactions are weak, es < 0.5e, 
a lubrication layer is formed and the slip length increases, 6 oc {rj/rjs) (^ + ^s)- (iii) If 
the coupling between solid and liquid is strong, however, the ratio rj/rjs decreases upon 
cooling and so does 6- If the ratio becomes sufficiently small (cf. Eq. (I33I1). as it does in 
the case e^ = 1 for a melt at an attractive surface, the Navier-slip condition fails. 

It is an open question by what condition the Navier-slip condition, Eq. ([2]), has to be 
replaced in case it fails. The schematic, two-layer model can account for the phenomenon 
that is observed in the molecular simulations but it introduces two additional parameters 
- the effective viscosity, ?7s, and the thickness. A, of the boundary layer. Moreover, the 
sharp interface between the boundary region and the bulk is only a crude representation 
of the rather gradual changes of the liquid structure and its dynamics in the vicinity of 
the surface. The stress balance between the viscous stress in the liquid and the friction 
stress at the solid substrate, which is the physical rational for Eq. ([2]), is also part of 
our schematic, two-layer model because we require the stress at the interface between 
boundary layer and bulk to be equal, cf. Eq. ( l27l) . Thus, one might speculate that it is 
the right hand side of Eq. ([2]) that relates the velocity gradient to the stress inside the 
fluid, which should be generalised in order to parameterise the changes of the fluid at 
the solid substrate. The failure of the simple, linear relation between stress and velocity 
gradient has often been observed for large shear rates [63l |64], where the polymer liquid 
exhibits shear-thinning. The schematic model, however, demonstrates that the shear 
rate and the concomitant non-Newtonian behaviour is not responsible for the failure 
and it suggests that the right hand side of Eq. ([2]) should be augmented by spatial 
derivatives of higher order. 

3-3- Soft, penetrable, brush-coated substrates 

The equilibrium and kinetic properties of the solid-liquid interface - wettability and 
hydrodynamic boundary condition - depend on the microscopic structure and dynamics 
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of the liquid at the sohd substrate. For instance, surface coatings - hke an oxide 
layer - can modify the wettability of a surface [THl |79] . Molecular simulations are well 
suited to investigate how one can tune these properties by controlling the microscopic 
interactions at the contact of solid and liquid. Since we utilise a coarse-grained polymer 
model [17], we focus the effect of a polymer brush coating. Polymer brushes offer a 
very stable and versatile strategy of tuning wettability, adhesion and friction [801 111]. 
Moreover, the effect of a polymer brush is rather universal because it much depends 
on the interdigitation between the brush and melt. The structure of the brush-melt 
interface is dictated by a balance of translational entropy, which favours mixing between 
brush and melt chains, and the loss of conformational entropy of the tethered chains, 
which stretch to allow the melt chains to come closer to the surface [81] |82| ^^. Both 
effects do not depend on the details of the interactions on the atomistic scale and can be 
described by coarse-grained polymer models. Additionally, the two contributions to the 
free energy can be controlled in the simulations as well as in experiments by changing 
the molecular weight of the free molecules of the liquid or the tethered chains of the 
brush and the grafting density of the brush. 

In the following we investigate a deceptively simple system comprised of polymer 
chains irreversibly tethered with one end to a solid substrate and a polymer melt of 
free polymer chains with the same molecular weight and interactions. In this way, the 
chains of the liquid and brush are of identical chemical nature. The grafting substrate 
is perfectly fiat and all friction between the substrate and the polymer melt stems from 
the interaction between the grafted chains of the brush and the liquid. 

Varying the grafting density gives rise to a rich wetting behaviour [831 EH ESI EH [IS] . 
Three different regimes can be distinguished: (i) At very low grafting densities, 
wettability is controlled by the interactions between the solid grafting substrate and the 
liquid. If the solid substrate is sufficiently attractive, the liquid will wet the substrate 
and this wetting transition is of first order [84]. Increasing the grafting density, the 
tethered chains provide additional attractive interactions for the liquid and the wetting 
transition occurs at a smaller value of the Hamaker constant. A, that characterise the 
strength of the long-range interaction between the solid substrate and the liquid, (ii) 
At intermediate grafting densities, the tethered chains provide enough attraction for 
the liquid to wet the brush even in the absence of additional attraction between the 
solid and the liquid. Long-range, second-order wetting transition can be observed as 
A changes its sign [84]. (iii) At high grafting densities, an interface between the brush 
and the melt gradually builds up and this interface becomes the narrower the larger the 
grafting density becomes. In this limit, the free chains of the melt cannot penetrate 
the brush and do not benefit from the attractive interactions between brush and melt. 
The brush-melt interface is characterised by a thermodynamic interface tension and the 
melt dewets from the brush although is comprised of identical constituents. This effect 
is denoted autophobicity [851 ES EZ] ■ 

The properties of the brush-melt interface have been studied by non-equilibrium 
molecular dynamics simulation of a Lennard- Jones, bead-spring model and dynamic 
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Figure 7. Flow of a polymer melt past a brusli-costed surface, (a) Poiseuille and 
Couette flows using dynamic SCMF simulations. The solid line indicate the total 
velocity profile, dashed lines depicts the velocity of the brush and the non-grafted 
melt, respectively, (b) The same flows observed by molecular dynamics simulations 
of the coarse-grained model described in the text, xp and xc indicate the position in 
the channel at which the fluid velocity is extrapolated to 0. The insets highlight the 
region of the brush-melt interface in which inversion of the flow is observed for the 
total velocity profile. Adapted from [13]. 



SCMF simulations of a soft, coarse-grained polymer model [HH [151 HI]- In figure [7] we 
present the velocity profiles obtained from non-equilibrium simulations of the Lennard- 
Jones, bead-spring model for Couette fiow and Poiseuille flow. First, we notice that the 
velocity proflle of Couette flow extrapolates to zero at a position, xq, which is well inside 
of the channel, Xq > 0. Thus, one cannot simply identify the hydrodynamic position 
with the location of the solid grafting substrate, Xb = 0, because this identiflcation would 
result in a negative value, 6 = —xq < 0, according to Eq. f l26|) . This negative value is 
inconsistent with the relation, 5 = 77/A, and the fact that both the shear viscosity, 77, and 
the friction coefficient, A must be positive. More generally, this observation highlight 
that both parameters of the Navier-slip condition, Xh and 6, have to be determined 
independently. 

For a brush-coated substrate, we expect that for low grafting densities the 
hydrodynamic position, Xb, is close to the solid substrate. At high grafting densities. 
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however, the hydrodynamic position, Xb, is dictated by the location of the interface 
between the brush and the melt. Quantitatively analysing the velocity profiles of Couette 
and Poiseuille flow, we again find that the parameters of the Navier-slip condition cannot 
be extracted because xp > xq (cf . figure [7]) for a wide range of grafting densities [15] . 
This effect is observed for both, the Lennard- Jones, bead-spring model and the soft, 
coarse-grained model of the dynamic SCMF simulations, which correspond to quite 
different situations. Qualitatively, this failure of the Navier-slip condition can also be 
rationalised by the schematic, two-layer model: The grafted chains of the polymer brush 
dangle into the polymer liquid and increase the effective viscosity, 77s ; of the liquid that 
penetrates the brush. Milchev and co-workers have recently quantified the increase of 
the near-substrate viscosity by defining an effective local viscosity as the ratio between 
the local stress and the gradient of the velocity [HHl l89] . 

Alternatively, the motion of the fluid through the polymer brush has been conceived 
as the flow through a porous medium, which is described by Brinkmann's equation 
[90l I9T] . Describing the polymer brush as a static porous medium, which only imparts 
additional friction onto the melt, one would expect that the tethered chains tilt in the 
flow and slightly fluctuate around their average, tilted positions. In the simulations, 
however, we observe strong, non-Gaussian fluctuations of the molecular orientation in 
shear [13] . Similar, strongly non-Gaussian distributions of orientations have also been 
predicted for isolated, tethered chains under shear [92 j . 

Only for extremely small grafting densities, Ug, we observe slippage because the 
solid substrate is assumed to be perfectly flat and in the limit cxg — )■ we observe plug 
flow with an infinite slip length, 5 [31]. In the opposite limit of very high grafting 
densities, the interface between brush and melt becomes very narrow and the brush 
resembles a very dense solid with a rather smooth surface. In this extreme limit, we 
observe a finite shp length, 5 [15] . 

While the velocity profiles at the centre of a brush-coated channel exhibit 
the expected linear and parabolic velocity profiles for Couette and Poiseuille flow, 
respectively, we observe a reversal of the flow direction in the vicinity of the grafting 
surface. This small reversal is observed both in the dynamic SCMF simulations, which 
use rather small Weissenberg numbers, Wi, and large invariant degrees of polymerisation, 
M, as well as in simulations of the Lennard- Jones, bead-spring model, which utilise larger 
Wi and much smaller M but capture the non-crossability of the flexible, linear chain 
molecules and the changes of the fluid structure in the vicinity of the solid substrate 

The reversal of the flow direction is observed at intermediate grafting densities 
on the brush side of the well-developed brush melt interface. In this region, the total 
velocity of the system is dominated by the velocity of the segments of the brush. Since 
the brush is irreversibly grafted to the solid substrate, which is at rest, the time average 
of the velocity of each brush segment in the steady state has to vanish. Nevertheless, 
the motion of the brush segment can couple to the rotational component of the shear 
flow at the brush-melt interface, which imparts a correlation between the tangential 
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Figure 8. /|| = 0.04e/iT The chains are comprised of A^ = 32 segments that interact 
via a repulsive Lennard- Jones potential with Vc = 2^/^cr. ksT/e = 1.68, vM ~ 6 and 
(JgR^o — l-SA/". The four panels depict a sequence of conformations extracted directly 
from a trajectory of the simulations, (top-left) Thermal fluctuations expose the chain 
to regions of higher velocity and (top-right) the convection in the flow stretches the 
conformation parallel to the substrate, (bottom-right) When a fluctuation brings the 
chain closer to the substrate, it is exposed to a smaller or even reversed flow and 
(bottom-left) its lateral extension collapses. 



velocity of a grafted bead, v^^^^^{x), and its distance, x, from the grafting surface. If a 
brush segment is located far away from the grafting surface, it is exposed to the flow 
and partially convected resulting in v^^'^^ix) > for large x. Since the time average 
of f m"'"'^'^ vanishes, this must be compensated by a negative v^^^^^ix) when the brush 
segment is located closer to the grafting substrate, i.e., small values of x. 

This cyclic motion of grafted polymers in shear flow has previously been predicted 
and observed for isolated, tethered molecules [93l HU [92] . In our simulations we observe 
this cyclic tumbling motion in a polymer brush [13]. All chains perform this tumbling 
motion albeit in an unsynchronised manner. Their collective behaviour results in the 
reversal of the total flow direction in the vicinity of the grafting substrate. 

We expect that this inversion of the velocity is characteristic for intermediate 
grafting densities. At very small grafting densities, the flow at the surface is dominated 
by the velocity of the solvent. The cyclic motion of the grafted chains will exist but 
it will not be sufficient to invert the total velocity, which is dominated by the flow of 
the solvent. Contrary, if the brush was strongly stretched, the width of the brush- melt 
interface would be narrow, the height fluctuations of brush segments would be small, and 
- in the Lennard- Jones model and experiments - topological constraints would become 
more important. These effects are expected to shift the inversion zone away from the 
grafting surface and tend to reduce the inversion of the flow. 



4. Polymer drops on surfaces 

The properties of the confining boundaries exert a pronounced infiuence on the fiow 
behaviour at the ultimate vicinity of the surface. This is particularly true for nano- 
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and inicro-fluidic devices or small drops because of the larger surface-to-volume ratio. 
Boundary effects are difficult to predict in the framework of a continuum theory because 
particle properties, details of the interactions between particles and correlations become 
important {e.g., layer effects of a dense liquid at a surface). These local properties can 
be accurately and efficiently investigated by molecular dynamics simulation. 

In the section, we will explore, which parameters have to be determined from 
a microscopic, particle-based model, in order to describe the behaviour of drops on 
surfaces with slip. On the one hand, such a parameter-passing between molecular 
simulation and a continuum description is necessary in order to extrapolate the results 
of the molecular dynamics simulation to experimentally relevant length and time scales, 
which are difficult to assess by molecular dynamics simulations. On the other hand, 
the comparison between molecular simulations and continuum description using the 
independently determined parameters demonstrates that the continuum model captures 
the underlying physics. 

4.1. Static properties on hard and soft, deformable substrates 

First, we will consider the equilibrium properties of polymer droplets on a hard, 
corrugated substrate at high temperatures, ksT/e = 1.2 (cf. section [3?T|) . In this case, 
the hydrodynamic position agrees with the location of the interface between solid and 
liquid, and the slip length can be controlled by the attraction, es, between solid and 
liquid. The macroscopic, equilibrium shape of a droplet is dictated by the volume, V 
of the liquid and its contact angle, G. The latter quantity can be estimated from the 
interface and surface tensions according to Young's equation. 

Directly measuring the contact angle of microscopically small drops in the 
simulations, however, one faces multiple difficulties [33]: (i) The liquid- vapour interface 
in the vicinity of the contact line is strongly distorted. In a cap-shaped drop geometry, 
these effects are particularly pronounced because the tension of the three-phase contact 
line between solid, liquid and vapour deforms the droplet. To leading order, the line 
tension. A, gives rise to a dependence of the contact angle, 9, on the droplet radius, R, 
which takes the form of Gretz's equation [95l [96] 



7lv cos e + 7ls - 7vs + p . ^ = (34) 

it smB 

(ii) The curvature of the liquid-vapour interface and the concomitant Laplace pressure 

shift the liquid-vapour equilibrium from its bulk coexistence value. Moreover, only in 

a certain range of system sizes equilibrium droplets of a specified radius, R, can be 

observed. |97]. (iii) Microscopically, the contact angle is not well defined but it always 

involves an extrapolation of the position of the liquid- vapour interface towards the solid 

substrate. For small droplet sizes such an extrapolation imparts a significant error onto 

the estimates of the contact angle, G. 

Therefore, we have decided to study cylindrical droplets as depicted in figure |9l 

which span the simulation cell in one direction. The length of the contact line is always 
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2Lz independent from the size of the droplet [981 EH US] • This largely reduces the effect 
of the line tension, although minor effects due to the interaction of the liquid-vapour 
interface with the solid substrate remain. Moreover, the size of the cylindrical droplet, 
R ~ n^/^ increases faster with the number of particles, n, than for a spherical drop, 
R ~ n^/^. By virtue of the low temperature, vapour density and liquid compressibility 
are small, and the shift of the pressure away from the bulk coexistence value does not 
give rise to a noticeable change of the liquid or vapour properties. The third caveat can 
be mitigated by avoiding to extract the contact line by the details of the liquid-vapour 
interface at the three-phase contact but rather fitting the shape of the entire drop to 
the macroscopic prediction. For instance, one can measure the radius, R, of the drop 
and the height xcm? of its centre of mass as a function of the particle number, n. For a 
cylindrical drop, R and xcm are related to the contact angle, 9, via [TB] 



xcM 4 sin^e 

-^ = 3 2e-sin29 '^"^^ ^^'^ 

Other quantities, like the volume, V, of the drop and its contact area, A can be expressed 
as 

V 26 -sin 20 

(36) 



/?2L, 2 



■'z 



^ =2 sine (37) 
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Fitting those relationships as a function of n, one can rather accurately determine 
the contact angle, 9, and the droplet radius, i?, from molecular simulations of 
microscopically small, cylindrical droplets. The translational symmetry along the z- 
axis allows to average profiles and increases the statistics of the simulation data. 

If the substrate is soft, then the forces that act at the three-phase contact line 
will also deform the substrate and lift up the three-phase contact line. This effect is 
well-known for liquid substrates and it is also notable for soft solids, like brush-coated 
substrates or polymer networks. The free energy penalty of a surface deformation of 
molten brushes has been considered within the strong stretching approximation. These 
calculations typically consider an incompressible brush and a sharp brush-melt interface. 
Under these conditions height fluctuations of the brush with long lateral wavelengths will 
be strongly suppressed because of the prohibitive costs of stretching of the lateral chain 
conformations |100[ 1101] 1102] . In our simulations, however, the polymer liquid exhibits a 
non-trivial equation of state with a finite compressibility. Therefore, a deformation of the 
surface of the brush can result in density fluctuations and even a spatially homogeneous 
displacement of the brush-melt interface is possible [15] . 

The deformability of the brush surface decreases with the grafting density, pg, 
because (i) the number of chains, which are stretched, increases and (ii) the density 
inside the brush increases and its compressibility decreases in turn. We have tried to 
study the deformability at the three-phase contact line of a polymer brush and a droplet 
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Figure 9. Snapshot of a cylindrical polymer drop on a polymer brush. Brush segments 
are coloured in grey, grafted ends are depicted in light green. Segments of the melt 
are coloured in red. The figure corresponds to a snapshot of the simulation box for a 
grafting density PgR^o = 9.5, chain length iV = 40, and ksT/e — 1.2 (with Tc = 2.5ct). 




Figure 10. Density contour plot of a A^ = 10 mers polymer liquid droplet wetting 
a deformable substrate. The reduced grafting density of the polymer brush substrate 
is Pgi?eo = 13.1 and the brush-melt compatibility is ebm = 0.6e. The dashed line 
represents the average free surface position of the unperturbed brush. The lifting-up 
of the contact line and the formation of a ridge is visible. 



that is comprised of identical chains. In order to achieve a contact angle of about 90°, the 
grafting density has to be large. Under these conditions the brush is hardly deformable 
and resembles a solid substrate. If the grafting density is low, in turn, the brush will be 
deformable but the polymer liquid will wet the brush or will make a very small contact 
angle. In this case, the forces at the three-phase contact line are nearly tangential and 
the lifting-up of the contact line is minuscule. 

Therefore, we have decided to make the molecules of the brush and the polymer 
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Figure 11. Dependence of the shape of a polymer drop {N = 10) on a brush {N = 40) 
as a function of their compatibihty, ebm- The left panel (a) presents the contact angle, 
Q, for different grafting densities, while the right panel (b) depicts the height, h^dgo 
of the ridge at three-phase contact. 



liquid slightly incompatible by reducing the energy parameter of the Lennard- Jones 
potential between segments of the brush and melt, 0.3e < ebm < 0.8e. Varying ebm and 
the grafting density, Pg, we can independently control the contact angle, 6, and the 
deformability of the polymer brush, and thereby identify conditions, where the brush is 
soft and the liquid-vapour interface is perpendicular to the substrate at the three-phase 
contact line. Figure [10] depicts the density contours of such an incompatible polymer 
drop comprised of chains with A^ = 10 on a slightly incompatible polymer brush. 

From these contour plots one can roughly estimate the contact angle, O, of the 
polymer drop on top of the brush. These data are presented in figure [TTk for a brush 
composed of chains of lengths, A'^b = 40, and a melt of short polymers, A^ = 10. As one 
increases the compatibility between the brush and the melt, ebm -^ e, the contact angle 
decreases. For the value, ebm ~ 0.6e, we achieve a contact angle of 90°. Interestingly, 
the data show that the contact angle is rather independent from the grafting density 
of the brush. We also find that G does not significantly vary with the length of the 
grafted chains (data not shown). For a brush with an intermediate grafting densities in 
a bad solvent, the structure of the brush surface is independent from pg and resembles 
a liquid-vapour interface. The independence of the contact angle therefore implies that 
the brush-melt interface does not significantly change with grafting density, Pg. Thus, 
it appears that the value of ebm ~ 0.6e is sufficient to dominate the tension between the 
brush and the melt and to outweigh the subtle entropic contributions that determine 
autophobicity and that do depend on grafting density and the ratio of chain lengths 
between the grafted and free chains [81]. 

The height, /iridge, of the ridge as a function of the compatibility, ebm, is presented 
in panel (b) of figure [Til It is only of the size of the segment diameter, a, and it slightly 
increases with reducing the grafting density because the brush becomes softer. The main 
deformation, however, occurs at the brush-melt interface and thus the dependence on pg 
again is rather weak. Although the contact angle decreases with ebm? and therefore the 
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tension of the liquid- vapour interface acts more tangentially to the surface, the height, 
bridge; of the ridge increases as we make brush and melt more compatible. Note also 
that on the microscopic scale of the range of interaction, re = 2.5a, (cf. figure [TOj) . the 
liquid-vapour interface bents around before the ridge, where it meets the substrate. 

4-2. Motion of droplets on hard substrates 

To illustrate that the microscopic details of the structure at the liquid solid interface 
can be accounted for in a hydrodynamic description by the contact angle and the slip 
length, we have studied the motion of droplets under the influence of an external body 
force [16]. Examples of such a situation are falhng drops on an inchned plane due to 
gravity |103[ I104[ I105[ 1106] . Much attention has been devoted to the onset of motion, 
the changes of the shape of the droplets, and their instabilities in response to strong 
body forces |107[ I105[ I108[ 1109(1110] . Comparably less attention has been devoted to the 
steady-state rolling and sliding of drops in external fields |11H 1106] , which is pertinent 
to nano-fiuidic devices. 

The steady-state tells a great deal about the dissipation mechanisms inside the 
moving droplets for the drop's velocity is dictated by the balance between the power 
gained through the motion in the external field and the dissipation. The dominating 
source of dissipation depends on the size of the droplet |112j and a systematic study 
of the size dependence is indispensable for identifying the different contributions in a 
simulation study. In the following we use the same model as in section 13.11 on two layers 
of an FCC Lennard- Jones solid. 

Similar dissipation mechanisms also act when a droplet spreads on a wettable 
substrate [IIHl [IH [USl HTH HIZl HHl Hm |120] but the study of the steady-state 
properties is computationally much more convenient because one can average properties 
along the trajectory instead of averaging over many realisations of the spreading process. 

The total dissipation rate, T ^^Qr^ is additively comprised of different contributions 
[3]. Energy is dissipated in the flow of the viscous fluid inside the drop. Assuming 
that the velocity field is predominantly tangential to the substrate, the lubrication 
approximation estimates 

This contribution scales like the volume, V , of the drop. In the last step, we have rather 
crudely estimated the typical velocity gradient, -^ inside the drop by its value, -^, at 
the substrate according to the Navier-slip condition Eq. ([2]). ^ is a factor of order unity 
that describes the details of the velocity field inside the drop. The latter approximation 
is restricted to small contact angles and large slip lengths, 5. 

An additional source of dissipation are processes at the three-phase contact line 
[1211 11221 1123] . which inter alia may arise from pinning at the contact line or 
evaporation/condensation of liquid at the three-phase contact. This contribution, 
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however, scales only like the length of the three-phase contact line. In the case of 
cylindrical droplets it provides a constant dissipation independent of the drop size and 
even for cap-shaped drops we expect this contribution to only influence the behaviour 
of the smallest droplets. If the substrate was soft and deformable {e.g., a brush- 
coated solid), then the lifting-up of the three-phase contact and the deformation of 
the soft solid in the course of the motion would constitute an additionally dissipation 
mechanism at the contact line [124[ 11251 1126[ I127i 1128] . This effect has been termed 
"viscoelastic braking" |1 25] and is an interesting subject for future simulation studies 
using the substrates characterised in the previous subsections. 

Another contribution may arise from a precursor film that spreads in front of the 
drop. Dissipation in such a precursor film is important for spreading of drops but it is 
not observed in the present simulation of the motion of drops under the influence of a 
body force. Therefore it is not considered further in the following. 

If slip occurs at the solid-liquid boundary, however, an additional important 
dissipation mechanism arises that stems from the friction as the liquid slips past the 
solid. Since the entire contact area. A, moves, this dissipation mechanism scales like 
A and outweighs the dissipation at the three-phase contact. In the ultimate vicinity of 
the corrugated substrate, on the microscopic scale, the liquid flow follows the substrate 
shape and the concomitant velocity gradients result in dissipation. On the one hand, this 
friction can be conceived as a viscous dissipation at the solid-liquid boundary. On the 
other hand, the physical mechanisms of friction at the solid-liquid boundary resemble 
the processes at the three-phase contact line. [121] Although the origin of friction at the 
solid-liquid boundary is complex and possibly diverse, its effect can be quantitatively 
described via the slip length, 6. The power dissipated by the friction is 

Ty^ = FmcV\\s with Ffric = 77— ^A (39) 

A ^ 

where the dynamic friction force, -Ffric, is given by the definition of the Navier-slip Eq. ([2]) 
[59] and f ys is the velocity of the liquid at the substrate. 

From dimensional considerations, Eq. ( 138|) and Eq. (!39|) imply that the dissipation 
in large drops is dominated by the viscous dissipation of the liquid flow inside the drop 
while for smaller droplets that friction at the substrate plays in important role. Since 
the slip length, 6, is the only length scale the ratio, R/S, separates large drops from 
small droplets. Using the approximation for the viscous dissipation in the volume of the 
drop, one can make further progress by balancing the power gained by the motion of 
the drop with velocity, u, in the external body force, Pl/|i, against the two dominating 
dissipation mechanisms 

p,V/„„ = 4^ + f,(^)V~„(^)\M + {n (40) 

where, in the last step, we have estimated the magnitude of the velocity of the drop, 
u, by the tangential velocity of the liquid inside the droplet at the height, xqm, of its 
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Figure 12. a) Velocity profile as a function of the height, x for different lateral 
positions, y, inside the droplet. The inset shows the sections, which have been used 
to calculate the velocity profile v^^{x,y). b) Velocity of the centre of mass per unit 
acceleration as a function of the number of monomers for es — 0.4 e (O = 130°, upper 
panel) and es — 0.8 e (6 « 29°, lower panel). The circles are the results of the 
molecular dynamics simulations and the dashed lines correspond to Eq. (j42[) using the 
independently determined values of rj and S. The error bars mark the uncertainties 
associated with these independently determined input parameters. The asymptotic 
limits for small and large droplets are indicated by the solid lines proportional to ^/n 
and n, respectively. Adapted from Ref. |16) 



centre of mass, u ~ "ynfa^cM; 
velocity, u is given by |§| 



f lis + -^xqm- This relation implies that the steady-state 



rju I SR for SA ^ V dominated by friction 
Pl/|| I R^ for 6A '^V dominated by viscous dissipation 



(41) 



A somewhat more quantitative description of the velocity field inside the drop can be 
obtained by assuming that the tangential velocity profile, v\\{x), only depends on the 
distance, x, from the substrate. In this case, the Stokes equation for a liquid under the 
influence of a body force predicts a parabolic velocity profile 



V\\[X) 



Pl/|| 



H - 



Ax 



Ax + 6H 



(42) 



where Ax = x — xi, denotes the distance from the hydrodynamic boundary position, and 
H is the height of the drop with respect to xi,. The coefficient in front of the quadratic 
term is dictated by the driving force, /y, and the remaining two coefficients have been 
obtained from the Navier-slip condition Eq. ([2]) and the observation that there is no 
tangential stress at the liquid- vapour interface, i.e., -q^ = for Ax = H. In figure 
[T2] we compare this phenomenological hydrodynamic description with our simulation 
data. In panel (a) we demonstrate that the tangential velocity profile, v\\{x), is indeed 

§ The velocity, u, and the body force, pl/||, can be quantified by the dimensionless capillary number, 
Ca= ury/7Lv, and the Bond number, Bo= PhfwR" /jly- Then Eq. (|41[) can be rewritten in the form 
riu/{pi^f\\R^) ^ Ca/Bo ^ 5/ R for small droplets while the ratio Ca/Bo approaches a constant for large 
drops. 
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Figure 13. Dewetting profiles for different equilibrium contact angles as indicated 
in the key. The length of the substrate is 2000(7 and its width 20cr. Each system 
is comprised of 300 000 particles. The successive profiles correspond to snapshots at 
5000 T time intervals. 



largely independent from y and that Eq. fj42|) provides a rather accurate description 
using the previously determined parameters, r] and b. Panel (b) depicts the dependence 
of the drops velocity, u ~ v\\{xqm), on the size of the cylindrical drop and compares the 
simulation results to Eq. ( |42|) without adjustable parameter. Quantitative agreement is 
found. The description also captures the counter-intuitive observation that the crossover 
from the friction-dominated behaviour of small droplets to the viscosity-dominated 
behaviour of large drop occurs at smaller sizes. Re ~ 5/(1 — cosO) for the larger 
contact angles, i.e., the crossover between surface-dominated and volume-dominated 
dissipation occurs at a larger volume, n^ ~ 160 000, for systems with a larger contact 
angle (es = 0.4e) than for flatter droplets ( n^ ~ 33 000, es = 0.8e). Although, at flxed 
volume, the flatter drop makes a larger area with the solid substrate, the velocity fleld 
inside a flat drop also strongly deviates from that of a solid-body rotation and therefore 
signiflcantly contributes to the dissipation. 



J^.3. Outlook on dewetting of polymer films 

In equilibrium a partially wetting liquid on a solid substrate has the shape of a drop with 
a contact angle, 6, determined by Young's equation, Eq. ([1]). However, the polymer 
liquid can be initially prepared as a metastable fllm of uniform thickness on the substrate 
{e.g., by spin-coating). In this case, two situations can be distinguished: (i) If the film 
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is thin enough, the fluctuations of the local film thickness in the film will spontaneously 
grow everywhere in the film. Through this spinodal dewetting mechanism |129] . the film 
will rapidly relax to its equilibrium configuration {i.e., a drop), (ii) A thick film will be 
metastable to fiuctuations, and the dewetting process can be initiated on relevant time 
scale only by defects or a rupture of the film. The simulation data in figure [13] illustrate 
this second dewetting scenario. In this case a hole is nucleated and its radius grows in 
time. 

Both, the shape of the rim that forms as the liquid retracts from the hole and the 
dewetting speed, by which the hole's radius grows, depend on the microscopic properties 
of the substrate ^M EHl [76l [77] , the initial film thickness [Ml [l32l M, [133] and the 
viscous properties of the liquid |134[ I135[ I136[ I137[ I138[ 1139] . The larger the equilibrium 
contact angle the faster the dewetting, and this effect is also observed in the data 
depicted in figure [131 Moreover, slippage of the liquid past the substrate determines 
the dewetting speed and rim profile |138[ [TH [751 US HZ]- It can be shown that for plug 
fiow the hole radius, R, increases as i? ~ t^/^, while for a stick boundary condition one 
finds R r^ t |130j . As for droplets, this will cause a transition from surface- dominated 
to volume-dominated regime during the growth of the rim. The presence or absence of 
slip also affects the shape of the rim; indeed the rim profile can decay towards the outer, 
unperturbed region with oscillatory patterns as observed in experiments and predicted 
by analytical models [76]. Thus, a great deal of information about slippage and the 
hydrodynamic boundary conditions can be extracted from the speed and the shape of 
the dewetting front observed in experiments. 

Molecular dynamics simulation permits to have a microscopic insight into the 
dewetting dynamics and compare to the predictions of analytical models using the 
independently determined parameters like viscosity, surface tension and slip length. Thus 
one can gauge the validity of the necessary approximations in the analytical description 
of the macroscopic behaviour and indicate oversimplification {e.g., neglect of thermal 
fiuctuations). 

Unfortunately the simulation of dewetting phenomena poses a huge computational 
challenge because they require the study of large systems with more than 10^ beads 
over long times in order to observe the transition from slip-dominated to viscosity- 
dominated dewetting behaviour. Using efficient simulation techniques such studies have 
just become feasible as illustrated in figure [131 Already the preparation of the initial 
starting configuration of a homogeneous, thick polymer film is computationally difficult. 
Knowing the single-chain statistics and by using optimised packing algorithms, one 
can overcome this obstacle [140j and prepare a thin, equilibrated polymer film that is 
comprised of 300 000 segments on a hard solid substrate with FCC structure. Periodic 
boundaries in the y and z directions of dimensions 2000 a and 20 a are applied and 
the unperturbed film thickness amount to /i = 12 a. By varying the strength of the 
substrate e^/e between 0.4 — 0.8 one controls the equilibrium contact angle and shp 
length, and consequently the dewetting speed. Once the film equilibrated, we remove 
polymers at the centre in a slab of length 10 a and start the simulation. After a fast 



Flow past hard and soft surfaces 32 

relaxation at the three-phase contact hne the film starts to dewet the substrate. We 
depict in figure [13] the successive steps of dewetting from which the time dependence of 
the hole growth and the shape of the rim can be analysed similar to the experimental 
studies [76]. This analysis is ongoing. 

One important difference between the experiment and the simulations, however, is 
already visible in figure [13] In the experimental system, the profiles of the rim of the 
hole are perfectly smooth and fine details of the rim shape can be analysed from a single 
snapshot. In contract to the spinodal dewetting process [141J, thermal fiuctuations are 
not of crucial importance once a hole has formed. The length scale in the simulation, 
however, is much smaller and on these microscopic length scales thermal fluctuations do 
matter. For a single snapshot, these thermal fluctuations of the liquid- vapour interface, 
ie. capillary waves, are difficult to distinguish from the shape of the rim predicted by 
analytical approaches without fluctuations. In particular, one would have to average 
the simulations over multiple realisations of the dewetting process in order to decide if 
the profile of the rim decays oscillatory towards the initial film thickness. 

5. Perspectives 

Molecular simulations can build a bridge between the microscopic fiow behaviour of 
particle-based models and a continuum description. While those multi-scale techniques 
have matured for equilibrium properties [43l 11421 11431 11441 1145] . the development of 
computational strategies that can describe the dynamics of polymer liquids across 
different scales is still in its infancy. On the one hand, no single computational technique 
can simultaneously address the length and time scales from the microscopic structure 
and dynamics at the boundary between the solid substrate and the liquid to the motion 
of liquid drops in nano- or micro-fluidic devices. On the other hand, the microscopic 
details of the interactions dictate the behaviour on the larger scales. For instance, the 
microscopic interactions between atoms determine the viscosity of a liquid. This shear 
viscosity, in turn, controls how much liquid will pass through a pipe or a microfiuidic 
channel in response to a pressure drop. In order to achieve this goal, one can use (i) a 
coherent multi-scale simulation scheme [1461 11471 11481 1149] , which explicitly couples the 
particle-based microscopic model and a hydrodynamic, continuum description or (ii) a 
parameter-passing technique. We have adopted the latter approach in this work. 

The underlying idea is that the complex, microscopic behaviour inside the liquid and 
at its boundaries can be captured by a small number of coarse-grained parameters that, 
in turn, enter the hydrodynamic description. These coarse-grained parameters depend 
in a complex manner on the microscopic interactions, and it is the basic assumption that 
the behaviour on large time and length scales depends on the microscopic properties only 
via these coarse-grained parameters. The bulk behaviour of a simple, one-component 
liquid is characterised by its density, pl, and its shear viscosity, rj. In order to 
describe the behaviour at its boundaries, we have used the contact angle, 9, (or 
the free energy cost of a unit area of surface) as well as the slip length, 6, and the 
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position, Xb of the hydro dynamic boundary. Young's equation, Eq. ([1]), and the Navier- 
shp condition, Eq. (|2]), and the fact that these parameters enter into the continuum 
description of hquids in terms of the Navier-Stokes equation or approximations thereof 
{e.g., lubrication approximation |150j ) have motivated this specific choice of parameters. 

We would like to emphasise that this specific choice of coarse-grained parameters, 
Pl; 0; Vj ^j ^^d x\^, that are necessary to transfer knowledge from the microscopic 
model onto the level of a continuum description is not unique. This list depends 
on the specific problem that one wants to describe. For instance, for strong shear, 
the non-Newtonian character of polymer liquids will become important, which may be 
captured by a dependence of the viscosity on the shear rate (shear thinning) or frequency 
(visco-elasticity) . Under strong non-equilibrium conditions the temperature may not be 
constant in the liquid and the heat flux has to be considered. In the simplest case this 
will add the thermal conductivity to the list of coarse-grained parameters. There may 
be also additional parameters required to capture the behaviour at the boundary. In 
order to describe the details of the liquid-vapour interface in the vicinity of the contact 
line, e.g., one needs the interface potential, g{h), which measures the free energy cost 
per unit area of placing the liquid-vapour interface a distance, h, away from the solid 
substrate. 

The role of molecular simulations in this scope is two-fold: (i) Measurement of 
these coarse-grained parameters for a specific microscopic model: Such studies are 
practically useful because they predict the behaviour of liquids. An understanding 
of the molecular basis of the coarse-grained parameters will contribute to tailoring the 
properties of boundaries between solid and liquids. This aim requires accurate and 
efficient computational techniques for calculating the bulk and boundary properties for 
molecular model of liquids via simulations. For instance, several simulation techniques 
for calculating the shear viscosity {e.g., Green-Kubo relation, Eq. fITT]) . and reverse 
non-equilibrium molecular simulations, Eq. ( ITOj) . |36] ) and extracting the parameters of 
the hydrodynamic boundary condition have been discussed, (ii) Testing the validity 
of the parameter-passing strategy: By a judicious selection of model systems that 
incorporate the relevant aspects of the physical behaviour on the molecular scale, one 
can quantitatively compare molecular simulations and continuum description on scales 
where both descriptions are valid. First, these studies provide information down to 
which length scale a continuum description is valid and how it breaks down. Second, 
they allow to explore to what extent the choice of coarse-grained parameters will be 
sufficient to capture the relevant physics i.e., ask the question: What are the relevant 
coarse-grained parameters that transfer the physical properties at the molecular level 
into the continuum model? This question can be only answered for specific systems. 

For the motion of liquid drops with a rather large slip length on a hard, corrugated 
substrate, our simulations suggest that the coarse-grained parameters, which we have 
listed above, yield an appropriate description of the microscopic physics because using 
the independently determined parameters in an approximate hydrodynamic description 
provides a quantitative description of the velocity profile inside the drop and its steady- 
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state velocity [16]. A detailed comparison with a numerically accurate solution of the 
hydrodynamic theory for the steady-state motion is certainly warranted to validate this 
result and to help extrapolating the results to larger, experimentally accessible drop 
sizes. 

In the example of a brush-coated surface [131 [iH] or a strongly attractive solid 
[H], our simulations suggest that the shp length, S, and hydrodynamic position, Xb, do 
not provide an appropriate description of the flow at the boundary. The simulations 
demonstrate that these two parameters are not able to simultaneously describe planar 
shear flow (Couette) and parabolic, pressure-driven flow (Poiseuille). This failure of the 
Navier slip condition can be "reproduced" by a schematic, two-layer model demonstrating 
that it is not related to a shear-rate dependence of the viscosity (shear thinning) or a 
time-dependent viscosity (visco-elasticity) . We speculate that it truly is a boundary 
effect and not related to bulk properties of the fluid {e.g., non- Newtonian character 
of the polymer liquid). It is rather tempting to attribute this failure to higher-order 
spatial derivatives in the relation between viscous stress at the substrate and stress at 
the boundary. 

Since the question about the appropriate coarse-grained parameters for fluid flow 
at boundaries can only be answered for specific systems there is a rich and wide-open 
spectrum of systems and processes to be investigated. In this manuscript we have 
restricted our attention to the most basic, coarse-grained models of polymer liquids [i.e., 
a Lennard- Jones bead-spring model and a bead-spring model with soft potentials) and 
the most basic flow phenomena (Couette and Poiseuille flow, and the steady-state motion 
of droplets in response to a body force). Of course, nano- and microfluidic systems offer 
a much greater variety of flow geometries {e.g., chemically or topologically structured 
surfaces), deal with more comphcated fluids {e.g., viscoelastic liquids), and often require 
to consider additional phenomena {e.g., spreading of drops or mixing in mult i- component 
liquids, evaporation effects). It is an open, challenging, and very important question to 
generalise the concept of parameter-passing to these more complex situations. 
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